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Abstract: In this article, we show that the double scaling limit correlation functions 
of a random matrix model when two cuts merge with degeneracy 2m (i.e. when 
y ~ x 2m for arbitrary values of the integer m) are the same as the determinantal 

formulae defined by conformal (2m, 1) models. Our approach follows the one 
developed by Bergere and Eynard in [2] and uses a Lax pair representation of the 
conformal (2m, 1) models (giving Painleve II integrable hierarchy) as suggested by 
Bleher and Eynard in [1]. In particular we define Baker- Akhiezer functions associated 

to the Lax pair to construct a kernel which is then used to compute determinantal 
formulae giving the correlation functions of the double scaling limit of a matrix model 

near the merging of two cuts. 
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1 Introduction 

It has been known for a long time that the study of random matrix models in different 
scaling limits gives rise to a great number of well-known integrable equations; both 
PDEs of solitonic type (KdV and, more generally, Gelfand-Dikii equations) and ODEs 
arising from isomonodromic systems (like Painleve equations). A key idea in these 
studies is the notion of spectral curves attached to algebraic equations P(x, y) = 0. 
The genus of the curve gives the number of intervals on which the eigenvalues of the 
matrices will accumulate when their size tends to infinity. It is well known that, in the 
generic case, the curve behaves like y ~ \[x — a near a branchpoint a (an extremity 
of an interval); the appropriate double scaling limit gives the celebrated Airy kernel in 
connection with the (1, 2) minimal model. But it may happen by taking a fine-tuned 
limit (see for instance [2j [5]), that the behavior near a branchpoint differs from the 
generic case and takes the form of y q ~ [x — a) p . In such a case, it is expected that the 
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double scaling limit is related to the conformal (p, q) minimal model. In their articles 
[2] and [3], the authors opened the way to rigorous mathematical proofs in order to 
establish that the correlation functions of the double-scaling limit of a matrix model 
are the same as the ones defined by determinantal formulae arising from (p, q) models. 
In their articles, they apply this method to all (2m + 1, 2) models, i.e. suitable limits 
of matrix models where the spectral curve behaves like y 2 ~ x 2m+1 near an endpoint. 
In this article, we will use the same method for the (2m, 1) case which corresponds 
to a point where two cuts are merging with a degeneracy 2m. For a generic merging, 
i.e. m = 1 it has been proven in jl] that the suitable double scaling limit of the 
matrix model is connected to the Painleve II equation. Some similar results have been 
established with the study of a suitable Riemann-Hilbert problem. For example the 
case of an even-quartic polynomial has been studied in [5J. It would be also interesting 
to derive similar results, for these kernels, as the ones proved in [10]. Here, using 
the approach of [2], we find, as expected, that the correlation functions of the double 
scaling limit of the merging of two cuts with degeneracy 2m are expressed through the 
Lax system of the Painleve II hierarchy (see [TT] and |17j). 

2 Double scaling limit in random matrices: the 
merging of two cuts 

2.1 Hermitian matrix models and equilibrium density 

It is well known in the literature that the study of the Hermitian matrix model with 
partition function: 

Z N = I exp(-iV r Ir {V{M)))dM (2.1) 
Jw N 

with an even polynomial potential 

2d 

V(x) = Y,Ux i (2.2) 

i=l 

can be reduced to an eigenvalue problem: A = {Xj, j = 1, N} for the matrix M with 
distribution: 



[ exp [2 V log\X j -X k \-NY2v{X j )]dX l ...dX N (2.3) 

J* N \ 1 s aTu^ at 7~f I 



l<j<k<N i=l 
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When N — » oo, the distribution of the eigenvalues on the line dv^ix) = pN{x)dx is 
defined (in the distribution theory sense) by the formula 



4>{x)du]si{x) ■ 

Z 



(2.4) 

J3 1 

same as the equilibrium density du eq (x) given by the limit of the empirical density: 



For any test function <f>(x) there is a weak limit di/^Cx) := lim du^(x) which is the 

TV— »oo 



1 N 

du eg (x) = Jim^ ^2$(x- Xj) (2.5) 

For details about the existence of the distribution limits, the equality between the 
equilibrium density dv eq [x) and dv^x) and the following characterizations we refer 
the reader to [9], [16]. Nowadays, many properties of the equilibrium density are 
known. For example, we know [18] that the equilibrium density is supported by a 
finite number of intervals [aj, bj], j — 1, . . . , q and that it is absolutely continuous with 
respect to the Lebesgue measure: 

I i 

dv eq (x) = p(x)dx = ——h(x)R 1 ^ 2 (x)dx, R(x)dx = TT ((x — aMx — bA) dx (2.6) 
2m - LJ - 

J=l 

where h(x) is a polynomial of degree 2d — q — 1 and R 1 ^ 2 (x) is to be understood as the 
value on the upper cut of the principal sheet of the complex- valued function R l l 2 (z) 
with cuts on J = U?=i[%)^']- m ^ ac ^ the equilibrium density dv eq [x) is completely 
defined by the knowledge of the extremities a/s and bj's and the unknown coefficients 
of the polynomial h(x). It has been proved [12J that such quantities are uniquely 
determined by the following set of equations: 

1. Connexion between h(z) and the potential V(2:)|^j 

Viz) = Pol (h(z)R 1/2 (z)) (2.7) 

z— >oo 

2. Residue constraint: 

Res (h(z)R 1/2 (z)) = -2 (2.8) 

z— >oo 

3. Integrals constraints: 

/ 3+1 h(z)R^ 2 (z)dz = , Vj G {1, . . . , q - 1} (2.9) 



3 Hcre and below we denote with Pol(A(z)) the polynomial part of the asymptotic expansion of 
A{z) at infinity 
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Note also that the relation between h(z) and V(z) (2.7) can be inverted by: 

Viz) 



h(z) 



(2.10) 



In theory, the previous set of equation is sufficient to determine the whole solution 
dv eq [x) but, practically, since the equations are highly non-linear, it becomes very 
hard to compute the unknown coefficients for two or more intervals or for potentials 
of degree higher than 4. Moreover, in some exceptional situations, the previous set of 
equations has multiple solutions. In such situations, the good solution is determined 
by a positivity condition: 



h(x) > 0, Wx e J = [J[aj,bj] 

3=1 



(2.11) 



When Vx G Uj=i] a i>M : h(x) > 0> the potential V(x) and the equilibrium measure 
dv eq [x) are called regular. Otherwise the equilibrium density is called singular and the 
corresponding potential is called critical, meaning that there is at least one point on J 
where the equilibrium measure vanishes. For a regular potential, the situation can be 
summarized with the following picture: 



Density of eigenvalues 
in the large N limit 



Near a generic endpoint, the density 
behaves like a square root 




Figure 1: Example of a typical eigenvalue density for a regular potential. The density 
is spread here in three intervals. 



2.2 Singular densities for the (2m, 1) case 

In order to study what happens at a singular density, one embeds the potential V(x) 
into a parametric family V(x, t) so that for some t = t c the problem is at the criti- 
cal potential: V(x,t c ) = V c (x). Then the interesting questions are to determine the 
asymptotics of the eigenvalues correlation functions when t — > t c . Indeed for t ^ t c the 
potential is regular and we are in the situation depicted in the figure above. Therefore 
one can define aj(t), bj(t) and h(x,t) determining completely the equilibrium density 
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for t ytz t c and study their limits when t — > t c . In matrix models, it is often interesting 
to study a modified version of the integral (2.1) by introducing a parameter T often 
referred as "the temperature": 



N 

exp(- — Tr (V{M)))dM (2.12) 



It turns out that T can be used as a parameter for the study of singular densities. In 
order to fit into our previous description, we need to introduce the following notation: 

V(x,T) = ^ (2.13) 

To the study of the (2m, 1) model, we place ourselves in the case such that at T = T c 
the potential V(x,T c ) = V c (x) becomes singular (see formula (2.29) below) and gives 
rise to a singular density defined by the following 2m singular density: 

p(x, T c ) = Pc {x) = ^-{x- be) 2m Vx 2 - b 2 = ^-h c (x)Vx 2 -b 2 (2.14) 

2l7l~ 2l7l 

with e G] — 1, 1[ representing the position of the singular point in the interval } — b,b[ 
supporting the distribution. For T ^ T c , we assume that the density is supported by 
two intervals ]ai(T), bi(T)[ and ]<22(T), b 2 (T)[ and define (note the normalization with 

£): 

p{x,T) = ^Hx, T)y/(x - ai {T)){x - 6x(T))(x - a 2 {T)){x - b 2 {T)) (2.15) 
Note that in order to recover our singular density at T = T c we must have: 

1. ai(T) -b 

T-5>T C 

2. b x be 

T-s>T c 

3. a 2 (T) be 

T-s>T c 

4. b 2 {T) -> b 

T->T C 

5. (x — be)h(x,T) — > Th c (x) 

The previous assumptions correspond to the merging to two cuts with degeneracy 
2m (order of the singularity). The most general case would be a singular point a with 
pl(x) ~ (x — a) p , (p, q) G N 2 , which is expected to correspond to the (p, q) minimal 
model (for q > 2 we are speaking about multi-matrix models). In our case the situation 
can be summarized with the following picture: 
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Density of eigenvalues 
in the large TV limit at 
critical temperature 




Figure 2: Example of a critical eigenvalue density; at the point be the density is 
singular and it behaves like (x — be) 2m . 

In [I], the authors studied the case m = 1 in details and conjectured some connec- 
tions with Painleve II hierarchy for higher m. 

2.3 Double scaling limits in matrix models 

In the study of matrix models, one is usually interested in the following functions called 

resolvents: 




(2.16) 



and in their cumulants, also known as correlation functions: 




Here, the brackets stand for the integration relatively to the probability measure 
Z~^dv^{x\ the Aj's are the eigenvalues of the matrices and the index c stands for 
the cumulants part defined as follows: 

<A X > = < A x > c 
< A X A 2 > = < A X A 2 > c + < A x > c < A 2 > c 
< AxA 2 A 3 > = < AxA 2 A 3 > c + < A X A 2 > c < A 3 > + < A X A 3 > c < A 2 > 
+ < A 2 A 3 > c < A 3 > + < A x >< A 2 >< A 3 > 
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n k 

<A 1 ...A n > = < Aj >= n <A u>c 

k=l h\Jl2"\Jh=J *=1 



The joint density correlation functions p n (xx, . . . ,x n ) can easily be deduced from the 
former correlation functions: densities are discontinuities of the resolvents and resol- 
vents are Stieljes transforms of densities. For example: 

w x (x) = [ ^Xdx' pxix) = — (wAx - iO) - wAx + zO)) (2.18) 

J X — X ' 2l1T 

Then we want to use a formal ^ power-series development which unfortunately is not 
necessarily well-defined for all matrix models. Indeed, if one is interested in convergent 
matrix models, then one must be sure that such a series expansion commutes with 
integrations. In general, this does not happen and solutions of the convergent matrix 
model differ from the solutions of the formal matrix model (where by definition the 
development is assumed to exist and to commute with integrations). The explanation 
of this phenomenon is simple: when we use a series expansion, it automatically ignores 
the exponentially small factors (one can think, for example, to exp(— x 2 ) which has 
at x = oo the same asymptotic expansion as the zero function). To sum up, formal 
matrix models are easier to handle, because by definition the formal expansion exists 
and we can perform formal operations on it; but the price to pay is that we only 
get a part of the convergent solutions (we miss the exponentially decreasing terms). 
It could appear disappointing to consider just formal matrix models, since they do 
not carry the whole convergent solutions (and thus leads only to a significative but 
incomplete part of the convergent solutions), but fortunately differences between formal 
and convergent matrix models have been well studied, and in [6] , [32] , the authors show 
how to reconstruct with theta functions the convergent solutions from the formal ones. 
Moreover recently, in the article [7], the authors showed that once we assume the 
existence of the asymptotic expansion for the resolvent and the two-points function 



(formulas (2.20) below for n = 1,2), then the same kind of expansion is guaranteed 
for all correlations functions and for the partition function as well. From now on, we 
will place ourselves in the case of formal matrix models, i.e. we assume that there 
automatically exists an expansion of type: 



2-2g 



lnZ ^ = J2(r) h ( 2 - 19 ) 

9=0 



and 

°° / N\ 2 ~ 29 ~ n 

w n (x 1 ,...,x n ) = ^2l — \ w^(xi,...,x n ) (2.20) 



9=0 
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The numbers f g are called symplectic or spectral invariants of the model (invariant 
relatively to symplectic transformations of the spectral curve). We refer the reader to 
the recent article [S] about the existence of such an expansion. The previous expansion 
can be understood as a large N expansion and therefore in the limit iV — > oo one 
expects that the leading value (g = 0) corresponds to the "real" large N limit of the 
model. In fact this intuition is correct and it has been proved that 



y(x) = iirpeq(x) = ^V'(x) - wf'(x). 



(2.21) 



This formula establish a direct link between the equilibrium density and the leading 
order of the first correlation function. The function y(x) (which is up to a trivial 
rescaling the equilibrium density) is often named the spectral curve of the problem. In 
our case, it satisfies: 

y 2 {x) = Polynomial(x) = — h 2 (x, T)(x — ai(T))(x — bi(T))(x — a 2 (T))(x — 02{T)) 

(2.22) 

This identity defines the algebraic spectral curve y 1 = P{x) where P is a polynomial. 
We remind the reader that Eynard and Orantin showed in (T3J that for any algebraic 
curve P(x, y) = we can associate some symplectic invariants f g and Wn\x\, . . . , x n ). 
Moreover, when the algebraic curve comes from a matrix model, these invariants are 



the same as the one we defined earlier in (2.16) and (2.17). 

In our case, the function y(x, T) = ^h 2 (x, T)(x — ai(T))(x — bi(T))(x — a 2 (T))(x — 
&2CO) depends on the temperature T and so are the corresponding invariants 
w n (xi, . . . ,x n ,T) and f g (T). When T — > T c it is known that V<7 > 1, f g — > 00 and 



that the correlation functions diverges. This is so because the expansion (2.19) reaches 



its radius of convergence in T. In order to recover finite quantities, one has to rescale 
properly the variables at T ~ T c . In our case we will prove that the good rescaling is 
given by: 

x i = be+(T-T c )&Z i (2.23) 

so that 

1 ;£,T) 



1/rescaled 



lim ij(fie + (T-T e )*..i 

T->T C 



T-T r . 



(2.24) 



and 



"'rescaled,?!^!' 



lim 

T^T C 



w<n 9 \be + (T — T c )^^i, ■ ■ ■ ,be + (T — T c )^£ n ,T) 



and 



rescaled,g 



(t - T c y 



lim (T-T c y^f g 



(2.25) 



(2.26) 
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are finite quantities and that the new w^ caled „(£i, . . . , £ n ) an d /rescaied, 9 are the spectral 
invariants of the rescaled curve descaled (£)• I n ^he g enera l context of matrix model, 
such a rescaling is called a double scaling limit since we have performed a double limit 
N — > oo and T — > T c so that N(T — T c ) 2m remains finite: 



InZ 



A? 



N 



2-2g 



N 



2-2g 



rescaled, g 



(2.27) 



12 It) f^J2\ T 

9=0 V 7 9=0 K ° 

From a geometric point of view, this double scaling limit corresponds to a local zoom 
in the region of the degenerate point be. The rate of the zoom depends on both the 
temperature T and the size of the matrices N so that N(T — T c ) remains finite. It can 
be illustrated with the following picture: 



Density of eigenvalues 
in the large N limit 
at critical temperature 




Figure 3: Example of a critical eigenvalue density near the critical temperature 

In the context of matrix models, double scaling limits are often very important 
because they are expected to give universal (independent of the potential) rescaled 
spectral curve and correlation functions related to (p, q) minimal models (and thus in 
our case the (2m, 1) minimal model). On the other hand, (p,q) minimal models are 
studied through string reductions of some well known integrable systems. In the rest of 
the paper, we will prove that, in the case of the merging of two cuts, the rescaled spec- 
tral curve corresponds to the spectral curve of the (2m, 1) minimal model. Then, using 
the method introduced by Bergere and Eynard in [2] , we prove that the rescaled corre- 
lation functions and the spectral invariants correspond to some "correlation" functions 
expressed with some determinantal formulae [3] for the (2m, 1) minimal model. 



2.4 The rescaled spectral curve in our 2m degenerate matrix 
model case 

In order to get the rescaled spectral curve, we need to perform a few consecutive steps. 
First we can express explicitly the corresponding critical potential corresponding to 
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p c (x) in 2.14 using 2.7 The computation is straightforward and uses only the general 



Taylor expansion of: 



Vi + x 1 + 2^1 1J n!fn-l)!2 2 «~ 1 



n=l 



(2.28) 



It gives: 




2m+l-j ■ 



n=l 



2m \ (-l)*(2n - 2)!e 2 ( m -") +1 --''6 
In + j - l) n\{n - 1)!2 2 "- 1 

(2.29) 



where E( 2m+ ^ - ) stands for the greatest integer lower or equal to 2m ~^ J . The critical 



temperature is given by: 



L2m+2 m + 1 



,2m-2n+2 



(2m)! 



(2.30) 



■J n!(2m - 2n + 2)!(ti - 1)!2 2 "" 1 

Then, we need to use some reformulations of conditions 2^ and 2.10 Indeed, it is 
known for a long time (a proof can be found in appendix A of [2] but the results were 
derived much before) and has been used intensively in [TJ that the set of equations 
2.8 and 2.10 leads to the following ordinary differential equations (sometimes called 



hodograph equations): 
d 



dT 
d 

dT 
d 

dT 
d 

dT 



ax(T) 
a 2 (T) 



(2.31) 



4( ai (T) - x Q (T)) 


h{ ai {T),T){ ai {T) 


-b x (T))(a x (T)-a 2 (T))(a x (T)- 
— x (T)) 


~ b 2 (T)) 


h{b x {T),T){b x {T) 


- a x (T))(b x (T) - a 2 (T))(b x (T) - 
4(a 2 (T) - x (T)) 




h(a 2 (T),T)(a 2 (T) 


-b 2 (T))(a 2 (T)-a x (T))(a 2 (T)- 
A(b 2 (T) - x (T)) 


-h(T)) 


h(b 2 (T),T)(b 2 (T) 


- a 2 (T))(b 2 (T) - a x (T))(b 2 (T) - 


bi(T)) 



where the point Xq(T) (b x (T) < Xq(T) < a 2 (T) ) is determined by: 

r 2(T) z - x (T) 

k{T) ai(T))(6 2 (T) -z)(z- a 2 (T)) 

This set of equations taken at T = T c for a x and b 2 gives: 

da x {T) 2 



dz = 



(2.32) 



dT \T=T C 

db 2 (T) 

dT \T=T C 



(1 _|_ tymjjlm+l 

2 

f-y _ e \2mJ,2m+l 
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(2.33) 



so that in a neighbourhood of T c : 



fl l( T ) T ~ Tc b n +e \2m b 2m+ 



-(T-T c ) + o(T-T c 



(2.34) 



bo(T) ~ b-\ 

21 ; t=t c (1 - e ) 2m b 2m+1 



(T - T c ) + o(T - T c 



(2.35) 

As mentioned earlier, we expect that the functions Oj(T) and bj(T) will be analytic 
functions of A = (T — T c ) u , where v is an exponent that we will determine later. 
Therefore we introduce the following notations: 



h(T) = be + aA + J2Kn^ n 

n=l 

oo 

a 2 {T) = k + 7 A + ^a 2 , n A n 

n=l 

oo 

x (T) = be + X A + J2 x n& n 

n=l 

oo 

h(z,T) = T{z-bt) 2m - 1 + P(z)A + Y,h n (z)A r 



(2.36) 



n=l 



where P(z) and h n (z) are polynomials of degree at most 2m — 2. In equations 2.31 



for a 2 {T) and b\{T) we see that the l.h.s. is of order (T — T^ v ~ x whereas the r.h.s. 
is of order (T — j'J-C 2 ™- 1 )^ Hence, to have compatible equations we must have, as 
announced in the previous subsection, that 



v 



1 

2m 



(2.37) 



The next step is purely technical and consists in proving that a = —7. Since it is only 
a technical point, we postpone this discussion in Appendix [Bj With the help of this 
relation we can now determine the rescaled spectral curve. 
First remember that for T 



T c , we have 2.14 



( z - be) 2m = h c (z) = Pol 



V'(z) 



(2.38) 
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For T ^ T c , reminding that V(z,T) = and that p(z,T) is defined with a factor ^ 



in 



2.15 (which will cancel the one of V(x,T)) we have: 
h{z,T) = Pol ' 1 <:) 



(2.39) 



We now use the fact that up to order A 2m_1 , both a% and b 2 are respectively equal to 
—b and b (2.34). Therefore we get: 



h(z,T) = Pol 



V(z) 



y/(z 2 -b 2 )(z-a 2 )(z-b 1 ) 
Then, from the definition of h c (x) we have that: 



O (A 2m ) 



(2.40) 



so that: 



[z - be) 



V'(z) 



2m 



Pol 



V'(z) 



T c z^oo \ yfz 2 — b 2 



Vz 2 - b 2 
Putting back this identity into 



2.40 



T c (z - be) 1 ™ + O 
and noticing that 



(2.41) 
only gives 



v /(« 2 -6 a )(«-o 2 )(«-6i) 

negative powers of z that will disappear when taking the polynomial part, we find 



that: 



h(z,T) 



T c Pol 

z— >oo 



T c Pol 

z— »oo 



(z - be) 



2m 



a/(^ - a 2 )(z - 6i) 

(z - k) 2 " 1 " 1 



+ O (A 2m ) 



(2.42) 



i , 2be-a 2 -b 1 , (be— a2)(be— &i) 
1 ~T z-be («-6e) 2 



+ O (A 2m ) 



We can now insert the Taylor series of the square-root: 

oo 



n=0 



(~l)"(2n)! B 

f n n2 2 2n 



(2.43) 



to get: 



Pol 

z— >oo 



(u - h^rn-i (, (-l) n (2n)\ { 2be-a 2 -b 1 (be - a 2 )(be - 

[ [ ' V ^ N) 2 2 2 " V (z-fre) 1 



+0 (A 2m ) 
(z - fre) 2 ™^ 

oo n , 
n=l fc=0 V 7 



•l) n (2n)! fn 



(26 e -a 2 -6 1 ) fc ((6e-a 2 )(6e-6 1 )) 



n—k 
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h(z,T) 



T 

± c 



(z - be) 2m - 1+k - 2n j + O (A 2m ) 

(2.44) 

Let's now introduce the following ensemble : 

I m = {( n , k) G (N* x N) / 2n - k < 2m - 1 and k < n} (2.45) 
Clearly I m is a finite set and we can rewrite the previous identity as: 

(z - be) 2m - x + 

(n,fc)e/ m v ' ; v 7 

(z - 6e) 2m - 1+fc - 2n ) + O (A 2m ) 

(2.46) 

We can now introduce the series expansion in A: 

2be -a 2 -b 1 = -(a + 7) A + O (A 2 ) 

and 

{ be-a 2 ){be-b l )= UA + f> 2 , n A") ( 7 A + f>, n AM = A 2 cry + O (A 3 ) 

\ n=2 / \ n=2 / 

Then we perform the rescaling 

z = be + A£ (2.47) 

We only need to take into account terms with degree strictly less than A 2m so that 
only a few terms remain: 



+0 (A 2m ) (2.48) 

so that: 



'-rescaled 



k ( \n—k c2m—l+k—2n 



(n,k)<=I, 



(2.49) 

Eventually we get the rescaled spectral curve by taking into account the trivial term 
R l l 2 ( Zl T) = y/(z - ai(T))(z - a 2 (T))(z - &i(T))(z - 6 2 (T)) with the rescaling 



2.47 



i^(6e + £A,A) = ^(be + £A - ai (A))(be + £A - a 2 (A))(6e + £A - &i(A))(6e + £A - 6 2 (A)) 
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bV^ly/(be + £A - a 2 (A))(6e + £A - 61(A)) + O ( A 2m ) 
i&AvT^V(f-a)(f- 7 ) + O (A 2 ) 



(2.50) 



so that: 
p(fe e + £A,A) 



V(e-a)(e-7) 



2;/( 



giving that: 
l/rescaled (0 



(2.52) 



2vr 

E i 7^§^Q c- 1 )^ +7) fc («7)- fc e 2m - 1+fc - 2 ")A 

-O (A 2m+1 ) (2.51) 



&vrv / T^ 2 v / (e-a)(7-0 
£2m-i + (-l) n (2n)! (n\ x . 

(n,k)el m 



(n\) 2 2 2n \k 



(-l) fe (a + 7) fe (a7) n ^ 



In the appendix [Bj we prove that a = —7 so that it eventually leads to: 



m—l 



2n) 



a 



2nt2m-l-2n 



n=l 



(2.53) 

We can even compute the precise value of 7. Indeed, using 2.45 to compute the 



leading term of the A-expansion of h(a,2,T) and putting it back into 2.31 (and using 



the fact that with the definition of x 2.32 we have X = when a + 7 = 0) we have 



a = -7 with 7 2m = a 2m 



Am 



[m 



|\222m+l 



m— 1 



(2n)! \ 6 2 (l-e 2 )(2m)! 



n=0 



(n!) 2 2 2r 



(2.54) 

In this case, introducing the new variable s by £ = 7s or equivalently 

2 = 6e + 7As (2.55) 
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we get: 



a 



"7 



1/rcscaled 



( S ) = 67T 7 2m v / r 3 ^v / ( l -s 2 ) 



m— 1 



,2m- 



n=l 



(2tQ! 
'n!) 2 2 2n ' 



2m-l-2n 



(2.56) 

Eventually 2.56 shows as expected that when performing a double scaling limit 
z = be+'yAs (with 7 a complex number given by 2.54 whose argument gives oscillations 



in the (Re(z) ,lm(z)) plane), we recover a universal curve. In the next section, we will 



see that this rescaled spectral curve 2.53 is exactly, (up to the trivial normalization 
factor by/1 — e 2 ) the spectral curve arising in the Lax pair representation of the Painleve 
II hierarchy with t m = 1, all other t/s (See next section for a definition) taken to zero 



and the identification uo(t) = 7 (coherently with 3.14). Before proceeding in the 



study of the Lax pair representation, we remind the reader that from general results of 



(<?) 

rescaled, n 



and 



Eynard and Orantin [T4] , the rescaled invariants and correlation functions w 
jrcscaied.g are automatically the symplectic invariants and correlation functions of the 
new rescaled spectral curve y re scaied(6) an d thus do automatically satisfied the famous 
loop equations 



3 Correlation functions and invariants arising in 
the Lax pair representation of the (2m, 1) mini- 
mal model 

In the previous section, we have found the rescaled spectral curve for a double scaling 
limit of a 2m degenerate merging of two cuts in matrix models. As conjectured in 
jl], we expect that this universal double scaling limit is connected to the Painleve 
II hierarchy. In order to prove this result, we will follow the approach [2] developed 
and successfully applied for the (2m + 1, 2) models. It consists in finding a natural 
spectral curve yhax{x) from a Lax pair representation of the hierarchy and check that 
it is equal to our rescaled curve defined in the previous section. Then from another 
work of Bergere and Eynard, [2] we can define from the Lax pair representation some 
new correlation functions w!f\xi, . . . ,x n ) and invariants F g by some determinantal 
formulae and a suitable kernel. In particular, they proved that these new functions 
do satisfy the same loop equations as our correlations functions. Eventually, with the 
study of the pole structure and we will end by proving that our new correlation 
functions wji (sci, • • • ,x n ) and invariants F g are identical to the rescaled ones defined 
in the previous section. 
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3.1 A Lax pair representation for the (2m, 1) minimal model 



In their paper [I], the authors claimed that a good Lax pair representation for the 
(2m, 1) minimal model should be given by a set of two 2x2 matrices TZ(x,t) and 
T>(x, t) satisfying the following Lax pair representation: 



Ndx 
Ndi 



V{x,t) = V{x,t)^(x,t) 
V(x,t) = K(x,t)y(x,t) 



where \l/(x, t) is a two by two matrix whose entries will be written as: 

#(x,t) = 



ip(x,t) 4>(x,t) 
if)(x,t) if)(x,t) 

and satisfies the normalization det^(x,t) = 1. 

The compatibility condition of the Lax pair is then: 



---V{x,t),K{x,t)--- 



(3.1) 



(3.2) 







(3.3) 



In order to specify completely the Lax pair, we need to impose some conditions about 
the shape of the matrices TZ(x,t) and V(x, t). In our case we will assume: 



K{x,t) 



x + u(t) 
-x + u(t) 



and 



with 



V(x,t) = J^tMx^) 



(3.4) 



(3.5) 



k=0 



V k (x,t) 



(3.6) 



-A k (x,t) xB k (x,t) + C k (x,t) 

xB k (x, t) - C k (x, t) A k (x, t) 

and A k , B k , C k are polynomials of x of degree respectively 2k — 2, 2k — 2, 2k. Note 
that in the literature one can find several different Lax pair corresponding to the 
same problem. Indeed any conjugation (change of basis) give equivalent matrices that 
describe the same problem but in different coordinates (see section [4]). In fact any 
equivalent Lax pair can be used since the quantities we will define later will be invariant 
from this choice. In order to have more compact notation, we will use the following 
convention: a dot will indicate a derivative relatively to t normalized by a coefficient 
1/N, namely: 

def 1 Of(x,t) 



ffaty 



N dt 



(3.7) 
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Putting back this specific shape of matrices into the compatibility equation gives the 
following recursion: 



A 
Ck+i 
Bk+i 
Ak+i 



0,50 = 0,^ = 1 
x 2 C k + R k (u) 
x 2 B k + R k (u) 
x 2 A k + -R(u) 



(3.8) 



where R k and R k are the modified Gelfand-Dikii polynomials given by the following 



recursion: 



u 

Rq(u) = u R Q {u) = — 

1 d 2 » 

R k +i(u) = uR k (u) - ^j^ R k{u) 
d R k (u) = u—R k (u) 



dt 



(3.9) 



It is then easy to see that the matrices 1Z(x, t) and T>(x, t) satisfy 3.1 if and only if 
u{t) satisfies the string equation (see details in jl].) 



^t k R k {u{t)) = -tu{t) 



k=0 



(3.10) 



which gives an explicit differential equation of order m satisfied by u(t) (since the 



polynomials R k can be explicitly computed from the recursion 3.9). In particular the 
case m — 1 gives Painleve II equation: 

d 2 u 
~dT 2 



-(t) = 2u 3 (t)+A(t + t )u(t) 



(3.11) 



where to is a free parameter that can be set to by a time-translation t = t + to- 
Remark : Seculiar equations 

As it is always the case for a linear differential equation, we can get a seculiar 
equation on tp(x,t) by combining the two components of the differential equation in 



t given by 3.1 In our case, we find that both ip(x,t) and (j)(x,t) are solution of the 

u(t)if)(x, t) 



seculiar equation: 



ip(x,t) 



(u 2 (t) -x 2 )ip{x,t) 



(3.12) 



x + u(t) 

which by a simple standard change of variable can be transformed into a Schrodinger- 
like equation. 
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3.2 Large N development 



From the fact that a dot derivative contributes with a factor jj, it is easy to see from 



the string equation 3.10 that u(t) admits a series development at large N: 



u 



(*) = E 



j=0 



N 2 



(3.13) 



Note : The fact that u(t) admits such a development in and not 4 comes from 
the fact that the modified Gelfand-Dikii polynomials -Rfc's are a sum of terms involving 
only even numbers of dots-derivatives (i.e. even power of jj). 



Putting back this expansion into the string equation 3.10 and looking at the power 
of iV° of the series gives us that Uo(t) must satisfy the following algebraic relation: 



(3.14) 

From that result, it is then easy to see that the matrices TZ(x,t) and T>(x,t) also 
admit a large N expansion: 



K{x,t) 
and 





—x + u (t) 



x + uo(t)\ l / o m(t) 







+ 



V(x,t) = J2 



N 2 \utit) 
V,(x,t) 



3=0 



TZj(x, t) 
N 2 i 



3=0 



No 



(3.15) 



(3.16) 



where the first matrix can be explicitly computed: 

V (x,t) = 



t + B + C 

t + B -C 



(3.17) 



with 



B n 



m j—1 
j=l k=0 



2(j-fc)' 



-i (2fc)l 
2 2k (k\) 2 



^o(t) 2fe+1 

«o(t) 2fe 



(3.18) 



fc=i 



It should also be possible to find equations defining recursively the next matrices 
TZj(x,t) and T>j(x,t) by looking at the next orders in the series expansion. But since 
we will have no use of such results we do not mention them here. 
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3.3 Spectral Curve attached to the Lax pair 



By definition, the spectral curve of a differential system like 3.1 is given by det(yld — 
T>o(x,t)) = 0, that is to say by the large N limit of the eigenvalues of the spectral 
problem (which we expect to give the large N limit of our matrix model). Note in 
particular that this definition is independent of a change of basis (conjugation by a 
matrix). From all the previous results, we can compute this two by two determinant 
and get: 



y 



(B + C + t) (B -Co-t) = Bl - (Co + t)' 



(3.19) 



Then, since we have xB = u (t) (C + t) from 3.14 and 3.18 a straightforward compu- 
tation gives the product as: 



VhSuX 



P(x,t) = (u (t) 2 -x 2 ) 



K j=l k=0 



2 2k (k\) 



-u (ty 



(3.20) 



In particular in the specific case where Vj < m : tj = 0, and t m — 1, we find that 
the spectral curve reduces to: 



Vj < m : tj = 0, t r 



Uo(t) 2 — X 2 



m— 1 

E 

fc=0 



X 



2(m-fc)-l 



(2k)\ 



2 2k (k\y 



u (t) 



2 k 



(3.21) 

As expected, with the identification Uo(t) = 7 we recover exactly the 
rescaled-spectral curve of our matrix model [2.53 



Note: In 3.20 we can see that the only simple zeros of P(x,t) are at x = ±uo(t). 
Moreover since the polynomial P(x, t) is obviously even and that there is no constant 
term in x in the sum, we get that P(x, t) has a double zero at x = and has double 
roots at some points ±Aj, i — 1, . . . , m — 1 

3.4 Asymptotics of the matrix ty(x,t) 

The next step in the method of [2] is to determine an asymptotic of the functions ip(x, t) 



and 4>(x } t). From the Schrodinger-like equation 3.12, we have a BKW expansion: 



«( I , 1 )= S ( Ilt )e*"(l + ^ + ^ + . 
Putting back into the seculiar equation gives the following result: 



1 / u (t) +x \ 2 eNf t ^ U 2 {tl) _ x2dt , A + ipifat) + 



V2 \u (t) 



x 



N 
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\/2 v«o(*)+^y \ AT 

i / u (t) -x y N ft y / U 2 (tl) _ X 2 dt , ( Mx,t) 



One can easily check that at dominant order in N the previous asymptotics gives 
det(\I/(x, £)) = 1 + (^). The next step is to transform the integration over t in the 
exponential as a integral over x by using the property of the spectral curve. Indeed, the 
spectral curve defines a Riemann surface which can be parametrized locally by x(z,t) 
and y(z,t) where z is a running point on the Riemann surface. Thus, the function y 
can be seen as both a function of (z,t) or (x,t). In order to avoid confusion here, we 
will write differently the function when it is seen as a function of (z, t) or as a function 
of (x, t) (we put a tilda for the function in (x, t) and keep y for the function of (z, t)): 

y(x, t) = y/P(x,t) = y{z{x, t), t) (3.23) 

Then, using standard chain rule derivation, one can compute: 

dy dx dy dx dy dx 

dz~di ~ ~diWz = ~didz~ ( ' ' 



From the expression of the spectral curve 3.20 (which gives explicitly y(x,t)) one can 
compute ||: 



dy xuo(d t u ) ^ x^-^^ky. 



2k 



ta s r* — -(^'<Ax 2 ^- i - k \2k)\2k 

+x(d t u )^ul-x^iJ2t j J2 22fc(fe \ )2 Mtf k 

x(d t uo) (2j)!(2j) j 

(3 ' 25) 



2 



To get the last identity, we have used the string equation 3.14 for Uo(t). Therefore by 
introducing the parametrization: 

z 2 = u (t) 2 - x 2 <=> x 2 = u {t) 2 - z 2 (3.26) 

one finds that: 



// A _ dx _ V^l-x 2 dx _ u (d t u ) 
X{Z ' t] ~ dz ~ x ' dt - x {6 - Z7) 



20 



so that eventually: 

dy dx dy dx dy dx 
dz dt dt dz dt dz 

The last identity can be rewritten as: 



x 



y/ u o~ 



-l 



(3.28) 



dy dx dy dx 
dt dz dz dt 

(3.29) 

and interpreted as the remaining of a non- commutative structure of [P, Q] = 4 in the 
limit N — > oo which in such situations often transform into a Poisson structure for 
y(z, t) -H- P and x(z, t) Q by simply replacing the commutator with a Lie bracket: 

{y{z,t),x(z,t)} = \ (3.30) 

With the help of this structure, we can get a reformulation of the integral: 

dy 1 



hence: 



and 



dt x'{z) 



d J x ydx 
dt 



\Jul{t') -x 2 dt' = J zdt 
Eventually we have the following large N developments: 

1 (u 



ydx 



ip(x,t) 
(j)(x,t) 
ip(x,t) 

4>(x,t) 



+ x 



Uo(t) - 


X 


( u {t) + 


\u (t) 




Uo(t) - 


X 


u {t) + 


X 


u (t) - 


X 


u Q (t) - 


X 



e -Nf*ydx 1 + 



N 

01 (X,t) 



Nj x ydx , + 



N 
N 



+ 



+ . . 



( ,v. r w, ( 1 , M*,t) , _ 



(3.31) 

(3.32) 
(3.33) 



(3.34) 



3.5 Kernels and correlation functions in the Lax pair formal- 
ism 

It was established in [3] that one can define a kernel K{x\,X2) and define from it 
(through determinantal formulae) some functions W n {x\, . . . ,x n ) that have nice prop- 
erties. In particular the authors showed in [3] that these functions do satisfy some loop 
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equations and thus are likely to correspond to our matrix model correlation functions. 
Following [3] we define the kernel by: 

K( Xl , x 2 ) = M x MM-tM*M (3.35) 

X\ — X 2 

Then we define the (connected) correlation functions by: 

W-i{x) = i)'{x)4>(x) - ${x)(j){x) (3.36) 



W n (x 1 ,.. .,x n ) 



S 



n,2 



(Xi - X 2 ) 2 

and eventually we define non-connected functions W nn - C by determinantal formulae: 



a=cycles i=l 



(3.37) 



W n ,„- C (xi, ...,x n ) = det(K(xi, xj)) 



(3.38) 



where the notation det means that the determinant is computed in the usual way as 
a sum over permutations o of products (—l) a Yl^ =1 K(x i ,x ai ), except for terms when 
% = cr(i) and when i = a(j) ,j = cr(i). In such cases, one must replace K(xi,Xi) by 
W\{xi) and K(xi,Xj)K(xj,xi) by —W 2 (xi,Xj). For additional details, we invite the 
reader to look at ([3]) 

As in our problem we will need the large N developments of these functions, we 
introduce the notations: 



K(x 1 ,x 2 ) 
W n (xi, ...,x n ) 



K (x 1 ,x 2 )e Nj ^ Hdx [l + ^2N- g K^( 



x u x 2 ) 



9=1 



J2N 2 - 29 ~ n w!f\x 1 ,...,x n ) 



9=0 



9=0 



(3.39) 



Then, we can insert all our previous results concerning the leading terms of the series 
expansion 3.34[3.37 and 3.39 It gives: 



K (xi,x 2 ) 



1 



i i i i- 

Uq + Xi \ 4 f u — X 2 \ 1 ( U - Xi \ 4 ( u + x 2 ^ 1 



2{X l - X 2 ) \ \U - X 1 J \Uq + X 2 



Uq + Xi 



Uq - X 2 



(3.40) 



W[°\x) = y(x) 



(3.41) 
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and 



4(xi - x 2 y 



-2 



(uq + xi)(u - x 2 ) I(uq -xi)(u + x 2 ) 



(u - xi)(u + x 2 ) y (u + Xi)(uq - x 2 ) 

(3.42) 

In order to get rid of the square-roots in the expressions above, it is better to 



introduce a proper parametrization of our spectral curve 3.20 Let us define: 




(3.43) 



In particular, under such a change of variables we obtain several useful identities: 

Uq — x .z — 1 



u + x 

Uq — X 
Uq — X 



Uq — X 2 

dx(z) 



dz 



z + 1 

u (z-l) 2 
2z 

Up(z + l) 2 
2z 

£(. + ix»-i> 

U (z 2 - 

2z 2 



1 



(3.44) 



Eventually we can rewrite W 2 ^ in terms of the new variable z: 

W 2 {0 \z u z 2 ) 



Aziz! 



(3.45) 



u^-l)(^-l)(^ 2 -l)^ 

Although these functions have some interesting features, they still depend on the 
choice of coordinates on the Riemann surface defined by the spectral curve. Therefore, 
we introduce similarly to [2] and [13] the corresponding differential forms: 



W£>(zi, ...,z n ) = W^\x{ Zl ), x{z n ))x\z x ) . . . x\z n ) + S n , 2 S g , , ^iMgg) 

\X\Z\) — x\z 2 )) 

(3.46) 

These differentials are symmetric rational functions of all their variables. Moreover as 



proved in the crucial theorem 3.2 these functions only have poles at z% = ±1 (except 
again (zi,z 2 ) which may have a pole at x(z\) = x(z 2 )). Eventually, a direct 
computation from |3.45 gives: 




(3.47) 
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3.6 Loop equations, determinantal formulae, pole structure 
and unicity 

The previous determinantal definitions may seem rather arbitrary, but as we mention 
before they have the interesting property (proved in [3]) to satisfy the following loop 
equations. 

Theorem 3.1 Loop equations satisfied by the determinantal functions: 

P-niX-i ■ • ■ ) Xn) ^^n+2,n— c\X j X, X\^ • • • j X n ) 

+ d W n (x,x 1 ,...,x j -i,x j+1 ,...,x n ) - W n (xx,. . . ,x n ) 



(3.48) 



JU >1 tAj JU 'J 

j=l J J 



is a polynomial of the variable x. The previous theorem is equivalently reformulated 
for the standard connected functions after projection on N 2 ~ 2g by the sets of equations 
(valid for every g > 0): 

h=0 ICJ 

A d wi 9 \x,J/{x J })-Wi 9 \x u ...,x n ) 

' J rp . nr> rp . 

\J >Xj -> JU JU i 

j — l J J 

(3.49) 

is a polynomial of the variable x. 

We emphasize again that loop equations are an essential step because it is well 
known in the matrix model world [IB] that the correlation functions introduced in our 
first section do satisfy these loop equations. Unfortunately, loop equations generally 
admit several solutions encoded essentially in the unknown coefficients of the polyno- 
mial P n . Therefore we need some additional results to get unicity. The first one deals 
with the pole structure: 

Theorem 3.2 Pole Structure : 

The functions z — > xjjk(z,t) are rational functions with poles only at z £ {±i, 0, oo}. 
The coefficients of these fractions depend on Uo(t) and its derivatives. Hence the de- 
terminantal correlation functions Wn^ are symmetric and rational functions in the 
variables z$ with poles only at z$ = ±1. 

Proof : The last part of the theorem is obvious from the definitions as soon as the 
results regarding the ipk(z,tys are established. This proof is presented in Appendix |A| 
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and is highly non-trivial. It uses the whole structure of integrability (i.e. the two dif- 
ferential equations 3.1 ) to eliminate other possible poles (at the other zeros of y-L^x)). 

With the knowledge of the pole structure of the Wn \ the fact that they satisfy the 
loop equations and the knowledge of we have a unicity theorem. In fact under 
these conditions we can identify our differentials Will's with the ones defined by the 
standard recursion relation introduced by Eynard and Orantin in [13]: 

Theorem 3.3 The differentials wi J ' satisfy the following recursion: 

VVSiOi, ■ ■ ■ , z n , z n+l ) = Res — — d z l '—j [W^iz, z,Z!,..., z n ) 

*-*±i 2u y(z){l - ^)(- - Zn+i) 

9 ' 



+EE<Vi^ / ) w iAi(^ J / J )] 



(3.50) 

3 ' 

where J is a short-writing for J = (z±, . . . , z n ) and ^ means that we exclude the 

h=oieJ 

terms (h, I) = (0, 0) and (h, I) = (g, J) in the sum. The notation z stands for the 
conjugate point of z near the poles where the residue is taken. In our case: z = - 

Note : It is worth noticing that in Eynard and Orantin's notation we have in our 
case (we omit the dependance in the t parameter): 

u(z) = y{z) 

_ z 

y( z ) = yi l / z ) = -y( z ) 

dE > { » ] = U*(I^ = 2( 2 -rt(^-l) (3 ' 51) 



so that: 

dE z (z n+1 ) z 2 



u(z) 2u y(z)(z - z n+1 )(l - z n+1 z) 2u y(z) (1 - **i)(± - z n+1 ) 



(3.52) 



Proof of 13.501 The unicity proof has been done in various article but for complete- 
ness we rederive it here with our notations. First of all Cauchy's theorem states that: 

W^i,...,4 + i)= Res dZ W^ 1 (z 1 ,...,z n+1 ) (3.53) 

We can move the integration contour to enclose all other poles, i.e. only ±1 in our 

case: 

W ( n 9 Uz U ---, z n + i) = Res - Wff !(*!,..., z ra+ i) 

z->±l z n+1 - z 
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z^±l Z n+ i - Z 



(3.54) 



We observe that the loop equations 3.49 can be rewritten in accordance with 3.49 
the following form by isolating the coefficients in the sum: 



in 



2W^(x)wii 1 (x l ,...,x n ,x) = EE<m(^ J )<n-|/|^ J / J ) 



h=0 ICJ 



d W ( n 9 \x, J/{xj}) - W^{Xj, J/{xj}) 



(9), 



dxj 



(3.55) 



Then putting this back into the residue computation, observing that the polynomial 
Pn\x; Xi, . . . ,x n ) does not contribute to the residue, and using the relation between 
x and z we are left with 13.501 



4 Lax pairs for the (2m, 1) minimal model and for 
the Painleve II hierarchy 

4.1 The (2m, 1) minimal model and the Flashka-Newell Lax 
pair 



As observed in [4] the string equation 3.10 is nothing but the m th member of the so- 



called Painleve II hierarchy. The Painleve II (PII) hierarchy, a collection of ODEs of 
order 2m, arises as a self-similar reduction of the mKdV hierarchy. In the papers [TT] 
and [IT] this relationship has been used to construct a Lax pair for the PII hierarchy 
starting from the relevant Lax pair for the modified KdV hierarchy. We call this PII 
Lax pair the Flashka-Newell Lax pair since the first member of the hierarchy was find, 
for the first time, in [15]. In this subsection we prove that, up to a linear transformation 
of the wave function and a rescaling of the variables, the Flashka-Newell Lax Pair is 
equivalent to the (2m, 1) minimal model Lax pair. In order to simplify notation we 
forget, in this section, the rescaling given by 1/N over the variables x and t. We begin 
with the case t± = = t 2 = . . ■ = t m -i- 

Proposition 4.1 Define ^ as a new wave function 

§ := jm 

with 



J :-- 



1 % 
i 1 



26 



and set t m i — > (4 m+1 /2) (all other parameters tj equal to 0). Then ^ satisfies the 
Flashka-Newell Lax pair as written in 07] /. 

Proof Since J is constant we observe that \P solve the Lax system 

= K(x,t)*(x,t) (4.1) 
with T>(x,t),TZ(x,t) obtained through conjugation with J; i.e. 

7t{x,t) = JKixrfJ- 1 = 

and 

f) (rr A — 4 m+1 TT) (rp +\ r-i — 4 m+1 ( —iC m (x,t) iA m (x,t) + xB m (x,i) 

u m {x,t)- 2 JU m {x,t)j - 2 y_ iAm ^ + xBmM iC m (x,t). 

These two matrices are exactly the ones appearing in (16a) and (16b) in [11J (modulo 
the identification u < — > w,x < — > X,t < — > z). For the matrix TZ this is self-evident. 
For T> m we just have to observe that it has the same shape as the matrix written in 
the right-hand side of (16b) (see eqs (14); in particular the polar part in (16b) is zero 
thanks to (14b)). On the other hand this conditions, plus compatibility condition, 
determines uniquely V m . 

Of course the result above is extended to the case in which all tj enter in T> just 
taking linear combinations of the matrices studied in the previous proposition. This 
has been done, for the Flashka-Newell pair, in [T7j (note, nevertheless, that there the 
spectral parameter is rotated; A — > —iX). Hence we have the following proposition. 

Proposition 4.2 Under a rescaling of all time variables tj — > ^-tj the (2m, 1)- 
minimal model Lax pair is equivalent to the Flashka-Newell Lax pair for the PII hier- 
archy. 



5 Conclusion and outlooks 

In section 1, we have established that the double scaling limit of a matrix model with 
a 2m-degenerate point can define a universal rescaled spectral curve y r escaied( a; )- 111 
section 1 we also reminded that the correlation functions and symplectic invariants 
Wn\xi, . . . , x n ) and f g can also be rescaled in a suitable way in order to give some new 
functions w^ caled n {x\, . . . ,x n ) and new symplectic invariants / rC scaicd,g corresponding 
respectively to the correlation functions and symplectic invariants of the rescaled curve 
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y resca i e d(x). Then, starting from a Lax pair of the Painleve II hierarchy and using the 
same method as [2] we have constructed a spectral curve yLax{ x ) which coincides with 
y T escaied( x ) f° r a natural choice of the flow parameters tj's. Finally, with the definition 
of a suitable kernel and determinantal formulae, we have defined in the same way as 
[2] some functions Wn^ having interesting properties (loop equations). Studying in 
details the pole structure and computing W2 (^1,^2), we have eventually shown that 
the function Wn^s are in fact exactly the correlation functions of the curve yh ax ( x )- 
Since the two spectral curves are the same, we have proved the statement: 

Theorem 5.1 The correlation functions (and spectral curve) of the double scaling limit 
of a 2m- degenerate merging of two cuts are the same as the functions Wn (and spec- 
tral curve) defined by determinantal formulae of the integrable Painleve II hierarchy's 
kernel. 

This result reinforces the links between double scaling limit in matrix models and 
integrable (p, q) minimal models. With this new result and the one of [2] , the two 
models are shown to be identical for (p = 2m, q = 1) and (p = 2m + 1, q = 2) (m G N*). 
However even if this identity is expected to hold for every (p, q), some complete proofs 
as the one presented here are still missing. Indeed, if our reasoning may seem easy 
to generalize for arbitrary value of p and q, the crucial theorem |3.49| establishing that 
the functions W^f 1 coming from determinantal formulae do satisfy the loop equations 
(proved in [3]) is only valid for q < 2 at the moment. Therefore a good approach to 
the generalization for arbitrary value of (p, q) could be to first extend this theorem for 
every (p, q) and then to use the method presented here to extend the result. 

Another approach could be to use this approach to study other integrable systems 
whose Lax pairs are known. Indeed, it is possible to perform the same method as the 
one presented here for any Lax pair. In particular, for every Lax pair, it would be in- 
teresting to analyse the associated spectral curve and the corresponding determinantal 
correlation functions. 



A Appendix: Pole structure for ipk(z, t) 



In order to use the unicity theorem 



3.49 



showing that the Wn ^s are the expected 
correlation functions, we need to precise the pole structure of the function ipk( x , t)'s and 
0fc(x,t)'s from which they are defined. In order to determine the functions ^(x, t)'s, 
one can insert the series expansion |3.34| into the seculiar equations. Since the case 
if)k(x,t) and 4>k(x,tys are similar (they satisfy the same seculiar equation), we will 
focus only on the -^(x, t)'s. The main issue of this appendix is that putting the large 
iV asymptotics of ij)(x,t) 3.34 into the seculiar equation a priori gives unwanted poles 
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at the zeros of y(x) for ipk{x,t) that we need to rule out. It is the purpose of this 
appendix to explain how this can be done. 



A.l Study of the differential equation in t 

From the fact that u(t) satisfies the string equation we remind the reader that we have 

EH 



3=1 



From this, it follows that ^ is: 

' at 

^ = 1 (A 2) 

dt P^uo) 1 ' ] 

Performing more derivations relatively to t can give the derivatives of Uo(t) to any order 
as a fraction whose denominator is always a power of Pq(uq). For example: 

d 2 u P»(uo) 



dt 2 (PoK)) 3 
dPuo iff'M (^M) 2 

(^oM) 4 + (^(uo)) 5 1 ' ] 

and so on. 

As a consequence, any power of any derivative of Uq remains a rational function of Uq 
with poles only at the roots of Pq(x). For example, expressions like + (^f ) 2 

will be rational functions of uq with poles only at the roots of Pq(x). 

Now, putting back the development of u(t) = Uo(t) + ^fW- + ^jW- + . . . into the full 



string equation |3.10| gives that any subleading order Uk can be expressed as a rational 
function of u with poles only at the roots of Pq(x). 

Eventually, inserting the shape of the function ip(x,t) into the seculiar equation 
and evaluating the order N~ k gives the following equation VA; > 2: 

d t %g d t g d t 2ip k _ 2 1 / d t u \ 1 ^-4 / d t u \ 

= ^-^ d '^- — + 2{— X ) k + 2^{— I ) t J' 



Oifj f 0,u \ , 1 / 0,11 \ ( <h 0,1'j 

V- 1 k-1 } 

fc-2 



¥ 2gh \u + xj t _! + 2 |j \u + xj^ \ 4,i gh + h 



i=0 



where we have written in short: 

if>o(x,t) 
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h(x, t) 
g(x,t) 



(A.5) 



u (t) + x\ 1/4 
u (t) - X 



and the notation f-^-), stands for the term in N k in the expansion of Note in 

V u+x J k u+x 

particular that these terms can be expressed as a fraction with poles at ito(t) + x = 
and at P' (%(£)) = (the last are independent of x). For example the first one is: 



dti)i{x,t) 



(d t u ) 2 x 2 (u - x)* Oo + x) 



4(w + 



+ 



[u - x 



-u 2 (t) 



(A.6) 



where remember that U2{t) can be expressed as a rational function of ito(t) whose poles 
are known to be only when uo(t) is at a root of Pq (and thus are independent of x). 
From this expression, it is clear that ipi(x,t) may only have x-dependent singularities 
at x = ±u and at x = oo. 

A. 2 Study of the differential equation in x 

The technic presented in the previous subsection can be carried out for the differential 



equation in x. Starting with the second equation of the Lax pair 3.1 
1 d 



-ty(x,t) 



-A(x,t) 



xB(x,t) + C(x,t) 



Ndx~ K ~'" J \xB(x,t) - C(x,t) A{x,t) 
we can derive another seculiar equation for both ijj(x,t) and 4>(x,t) 



V(x,t) (A.7) 







1 d 2 

iVW^'^ A' 2 



1 fd x {xB + C)\ d 



xB + C 
d x (xB + C) 



+ N[ d * A ~ A xB + C 



ip(x,t) -y 2 (x,t)^(x,t) 



(A.8) 



where we have used that: 



det(^) = 1 y 2 (x, t) = A(x, t) 2 + x 2 B(x, t) 2 - C(x, t) 2 (A.9) 
Note in particular in the last identity that the r.h.s. should have a large development 



whereas the l.h.s. y(x) given by 3.20 does not. Therefore, the l.h.s. must have vanishing 
subleading orders in -j^, Vk > 0. 



Moreover, reformulating 3.17 gives 



A 

(xB + C) 



y(x,t) 



Uq + X 

u — X 



(A.10) 
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where the subscript stands for the first order in the large N expansion. Indeed, it 
comes from the fact that: 



y 2 (x,t) 

(xB + C) = 

(xB -C) = P 2 (x,t) 

and eventually: 



{xB + C% {xB - C% = P(x, t) = Pt(x, t)P 2 {x, t) 



=1 fc=0 



2 2fe (A;!) 2 



(A.ll) 



v j=l fe=o 



2 2fc (£;!) 



P 2 {x,t) = Px{x,t) 



Uq — X 
U + X 



With [A7T0] it is easy to see that: 

'd x (xB + C)\ _ r7,// 



xB + C 



u 



y Uq — x 2 



(A.12) 



(A.13) 



which will be crucial for the coherence of the computation. Indeed, putting the large 
N expansion of if)(x, t): 

= 9 (M)e™<-> (l + ^ + + .. 



into 



A. 8 and comparing the first orders in gives: 








g(x, t)y 2 (x, t) - g(x, t)y 2 (x, t) 

d x (g(x, t)y(x, t)) + y{x, t)d x g(x, t) + g(x, t)y(x, t) 



d x {xB + C) 
xB + C 



(A. 14) 



The second equation with the help of A.13 determines g(x, t) coherently with A. 5 that 
is to say: 

'u (t) + x\ 1/4 



g(x,t) 



u (t) - X 



Note now that VA; > 0, the function only has singularities at the sineul 



ar- 



ities of 



i 



xBq+Cq 



according to the standard rules of Taylor series for a fraction. The next 



order, gives us the function ipi^x^t) (with the notation that a subscript k defines 
the term in N~ k in the expansion at large N): 



d x zg d x g ( d x (xB + C) 



2gy 2gy \ xB + C 



1 f d x (xB + C) 

2 V xB + C 



d r A - A 



d x (xB + C) 
xB + C 
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(A.15) 



From the definition of g(x,t), it is easy to compute: 



9x9 

9 

d x 2g 

9 



w 



2 Uq — x 2 



UqX 



+ 



dr. 



x- 



4K 



X 



2\2 



(A.16) 



and thus to see that d x tp\(x,t) is a function of x that may only have singularities at 
x = ±it , at a; = oo and at the others zeros of y(x) = 0. (it is so because { ^^b+c^ 
have the same singularities as xB ^ +c which by A. 10 are only at x — ±Uq, x = oo and 
at the zeros of y(x)). 

It is then possible to extend this result for higher terms in the large N expansion. 
The power gives: 

d x ig 



d x ip, 



k-i 



- — 7/> fc _2 - —d x 4> k _ 2 - -d x 2^ k _ 2 

2gy gy y 

'd x (xB + C)\ () r g v - 

k-i-i 2gy^ 



fc-2 



d x {xB + C) 
xB + C 



2 

i=0 v ' K-L-t fa i=Q 

+ 2y^{ xB + C ) k _J' A 2y^{ a ' A A xB + C J, 



k-2-i 

dJxB + C)' 



(A.17) 



where we have defined i/)q = 1. The precise form of the relation is mostly irrelevant, 
but the main fact is that if all the ipi{x, t) with i < k are assumed to have singularities 
only at x = ±wo, x = oo and at the other zeros of y(x) = 0, then the same is true for 
d x ipk by a simple recursion. 



A. 3 Pole structure of ipk(x,t) 



With the help of A. 4 and A.17 we are now able to prove that the only singularities of 
x !->■ i/>k( x i t) are at x — ±u and at x = oo. 

From A. 4 we have shown that d t ipk{x,t) can only have singularities at x = ±w (t), 
at x = oo and when uo(t) is at a root of Pq. But from A.17 we have shown that 
d x ipk(x,t) can only have singularities at x — ±u (t), at x = oo and at the other zeros 

"u (t) 2h 



m j—1 

of y(x) = given by x = Xi(t) solution of J2^jYl ~ 

j=l k=0 



2(j-fc)-l 



(2k)\ 



in 



3.20 



But 



2 2fc (fc!) 2 

these poles are incompatible with the former result. Indeed if ipk(x,t) had a pole at 
x = Xi(t), then dtipk{x,t) would also have a pole at x — \i(t), but we have shown that 
the only x-dependent singularities of d t ipk(x,t) are at x ± uo(t) or x = oo giving rise 
to a contradiction. Therefore: x h- >■ ip k (x,t) has only singularities at x = ±w 
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(square-root poles) and x = oo (poles) and in particular has no pole at the 
other zeros of y(x) = 0. This result is highly non trivial because we need to combine 
the two differential equations (i.e. the whole integrable structure) to get it. Hence, 
the structure of integrability seems to play an important underlying role in the pole 
structure and we can hope that such a result could extend to every integrable system. 

A. 4 Pole structure in the z variable 

In order to have only poles (and not square root singularities), we want to shift the 
former result to the z variable defined by: 

9 Un — X 1 — Z 2 

z 2 = — ^x = u - 2 (A.18) 



Uq + X 1 + Z 2 



Note that we have the identities: 



dx \ 1 z 2 

— = {d t u )— — 2 
ot 1 + z 2 

dx 4u z 

~d~z = ~{l + z 2 ) 2 
2u 



Uq + X 
U — X 



1 + z 2 
-2z 2 

l + z 2 i 

(—U ) 2 



l 

Z2 

Az 2 u 2 fl-z 2 ^ 

vM = (TT# p °((iT^ , " il) 

d*9, (l + z 2 ) 2 



[z,t) 



dx*9, x (l + z 2 ) 2 (l + z 2 ) 4 , k , 

— = ,22 + KA 2 4 A -19 

Note also that every polynomial in x will give a polynomial in frp-, that is to say 
a rational function in z with poles at z 2 + 1 = 0. 
The rules for derivation give that: 

d t Mz,t) = fyfcOM) + (A.20) 

d z ip k {z,t) = —d x il)i{x,t) (A.21) 



where all these terms are already known from the previous sections. If one uses A. 19 
and the remark that a polynomial in x will give a rational function in z with poles at 
z 2 + 1 = (and remember that functions A,B,C are polynomials in x), one can see 
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that the singularities of ipk(x,t) at x — ±u (square-root type) and at x = oo (poles), 
will transform into poles at z = x = —u Q ), z = oo x = u ) and z = ±i 
x = oo). 

Hence we have the final result: VA; > 0: the functions z — > ipk{z, t) are rational 
functions with poles only at z G {±z, 0, oo}. The coefficients of these fractions 
depend on uo(t) and its derivatives. 



B Appendix: Discussion about a = —7 

When computing the rescaled spectral curve in the matrix model double scaling limit, 



we need to find a relationship between a and 7 that are given by 2.36 

00 

6j(T) = be + a A + ^ b 1>n A n 



n=l 

00 



a 2 (T) = 6e + 7 A + ^a 2 , n A r 



(B.l) 



n=l 



where we remind that A = (T — T c )a^. A first argument in favour of the fact that 
a = —7 is the case when e = 0. Indeed, in such a case, the situation is fully symmetric 
around the singular point 0. Therefore, one expects the two endpoints b\(T) and a 2 {T) 
to be symmetric around x = for every value of T around T c . In such a case the identity 
VT ~ T c : ci2 (T) = — bi(T) gives a = —7. When e ^ 0, we can carry out a similar 
reasoning at first orders in A. Indeed, if we center the origin at be, then as we observed 
it several times, the endpoints ai and 62 can be considered to be respectively —b and b 
up to order A 6 . Therefore in the function R^{x) they only add a multiplicative trivial 
factor depending on e (y/l — e 2 to be precise) which will not change the symmetry 
around be of the endpoints a 2 and 61 at first orders in A. 

Eventually, another more explicit approach is to put the developments |2.36| into 



all the equations 2.10 2.8 |2.9 and 2.32 determining h(z,T), xq(T) and the endpoints 
a\(T) , bi(T) , a 2 (T) and b 2 (T). Doing so leads to an algebraic equation of degree 2m 
connecting a and 7: 

Q(«, 7 )=0 (B.2) 

with Q a symmetric, homogeneous polynomial of degree 2m. Unfortunately the system 
does not admit a unique solution as soon as m > 1. Indeed, although the solution a = 
—7 is always there, when m > 1 there are also other possibilities such as a = A7, A G C 
and 7 satisfying an equation of degree 2m with complex coefficients. Though it might 
appear surprising that the set of equations may have several distinct solutions (thus 
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giving several eigenvalues density), one must remember that they are some additional 
constraints for the solution. Indeed, if one wants to have a density distribution, it 
means that all quantities involved must at least be real and positive. Therefore only 
the solution a = —7 is possible. 

Note : In fact a and 7 are not necessarily well defined. Indeed, they are only 
defined up to a multiplicative (2m) th root of unity since the equation defining them 
is homogeneous of degree 2m. This is because the notion of A = (T — T c ) 2^ is also 
ambiguous, whereas A 2m , a 2m and 7 2m are well-defined quantities, (which explain why 
the development in a\(T) and 62 CO is well defined). Indeed, if one changes: 



2 i n 7T 
2m 



Vn G {1, . . . , 2m — 1} : A — > A = Ae 2m , a — )■ a = ae 2m and 7 — > 7 = 7e~ 

(B.3) 



then 2.36 remains unchanged. With the change £ — > £ = £e 2™ , the rescaled spectral 



curve remains unchanged. 
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